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Qj \ We consider vacuum tunnehng of a new kind where the false vacua are not 

>-Ch ■ translationally invariant, but have topological defects that break some of their 

translational symmetries. In the particular case where the topological defects 

'k>i I are cosmic strings, we show the existence of an 0(2) x 0(2) symmetric bounce 

^ • configuration that allows a semi-classical estimate of the rate of cosmic string 

induced tunneling. A method of reduction is then suggested for simplifying 
the computation of the bounce action. Some phenomenological applications 
are described. 
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1 Introduction 

First order phase transitions in field theories may sometimes take place through the nucle- 
ation of bubbles seeded by massive excitations [Q] . Of the several distinct processes that 
can lead to the seeding, a particularly interesting one is due to the instability of topological 
defects (like cosmic strings and monopoles) with respect to spontaneous dissociation. The 
phenomenon is best understood by an example. Consider the SU{5) symmetric GUT the- 
ory. The sequence of phase transitions SU{5) -^ SU{A) x f/(l) -^ SU{3) x SU{2) x f/(l) 
occurs as the universe cools from temperatures of 10^^ GeV. At the end of the first phase 
transition SU{4) x U{1) monopoles are formed. These monopoles may have a structure 
such that at their cores the vacuum expectation values (VEV) of the scalar fields re- 
sponsible for the second phase transition lie near the minimum of the effective potential 
of the SU{3) X SU{2) x f/(l) phase. As the temperature drops to the point where the 
S'f/(4) X U{\) phase is metastable compared to the S't/(3) x SU{2) x U{1) phase, the 
monopoles become unstable if it is energetically favourable for their cores to expand. In 
other words the monopoles themselves become bubbles of true vacuum inside a sea of false 
vacuum. It was shown in ref. that this phenomenon may precede the nucleation of 
the usual bubbles in the false vacuum and the phase transition may take place essentially 
through bubbles seeded by monopoles. 

An analogous situation is the case where the topological defects remain classically 
stable as the supercoohng proceeds but may decay by quantum tunneling. The result is 
the same seeding effect discussed above. However seeded bubbles compete with unseeded 
bubbles on the same footing; both are quantum (or thermal) processes and suppressed 
due to the usual barrier penetration (or Boltzman) factors and one must compute the rate 
of nucleation of the two kinds of bubbles to decide whether the seeding has an appreciable 
effect or not. 

The phenomenon is interesting in theories where the false vacuum is so long lived that 
the phase transition never takes place and the universe remains trapped in the supercooled 
phase forever. In the minimal supersymmetric model (MSSM), for instance, it may be 
possible that the true vacuum is charge and color breaking (CCB) and the existing, 
viable vacuum is only a false vacuum. Recently it has been shown that such viable false 
vacua can exist in other low energy supersymmetry breaking models [^] that also permit 
the existence of U{1) cosmic strings. In all these cases, the decay rate of the viable false 
vacuum provides interesting constraints on the model. However, the existence of topo- 



logical defects in the false vacua of these theories is also quite generic and independent 
constraints can be obtained by considering the nucleation rate of seeded and unseeded 
bubbles. While the tunneling associated with unseeded bubbles is well understood semi- 
classically and the rate can be found by computing the action of 0(4) invariant bounce 
configurations of the underlying Euclidean field theory, a similar semiclassical treatment 
is needed in the case of the seeded bubbles for a useful estimate of the corresponding 
tunneling rates. 

In this paper we treat a local region of space containing a topological defect as a transla- 
tionally non-invariant (non-homogeneous) false vacuum which, classically, corresponds to 
a well defined local minimum of the energy functional of the theory. We show that there 
are new bounce configurations (or instantons) associated with these non-homogeneous 
vacua that permit a semiclassical understanding of the tunneling seeded by the defect. In 
particular there is a bounce with a large symmetry that appears as the global minimum 
of the action in an open space of field configurations. We then devise simple numerical 
methods to compute a lower bound on the corresponding semiclassical tunneling rates. 
Although the new bounce configurations are more complicated than the usual bounces, 
we use a method of "reduction" to simplify the computation to the point where it is not 
much harder than the 0(4) symmetric computations for the usual bounces 0. 

We restrict ourselves to the zero temperature tunneling, since that is the most inter- 
esting case for theories where the false vacuum is viable and required to be very long 
lived. However, generalizing to the finite temperature case is straightforward and the 
computation of the free energy of the critical bubble is analogous to the computation 
of the bounce action in a 3 dimensional Euclidean field theory. Our treatment is also 
applicable to all types of topological defects except global texture (see [Q and references 
therein for a discussion of different types of defects), but for brevity we limit ourselves 
to the case of cosmic strings only. In this case the effect of the tunneling is the random 
appearence of bead-like bubbles along the length of the cosmic string, which captures the 
essential feature of tunneling seeded by extended defects. Note that cosmic strings also 
happen to be the most interesting type of defects from the point of view of cosmology (||^ 
and references therein). 

The paper is organized as follows. In section 2 we prove the existence of the new 
type of bounce solutions and show that they appear as global minima of the action in 
an open space of field configurations. In section 3 we indicate some phenomenological 
applications and effects of the seeded tunneling phenomenon. Finally in section 4 we 



suggest a method of "reduction" that makes the estimation of the bounce action feasible 
even in comphcated field theory models. 

2 The Bounce as a Global Minimum of Action 

The semiclassical vacuum tunneling rate was computed in ref. by a saddle point 
evaluation of the Euclidean path integral. In the case of a theory with a single real scalar 
field with a false vacuum at = 0/, the transition probability per unit volume per unit 
time was found to be: 
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where Se[4>] is the Euclidean action for the so called 'bounce' configuration and d^ = 
dfj_dfj_. The prime on the determinant indicates that the zero eigenvalues are omitted. 

We will use the Abelian Higgs model with a single complex scalar $ to make our 
arguments concrete. To develop the formalism we start with the Euclidean action: 



Sr = d'^x 
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where D^^ = {d^ - iA^)^ and F^^ = d^A^ - d^A^. 

Let us suppose that the potential f/($) has a global minimum at $( and a local 
minimum at $j 7^ 0. The false vacuum at $/ may permit the existence of cosmic strings. 
A false vacuum with a single cosmic string of unit winding number is the lowest energy 
state with the following boundary condition for the VEV of the scalar field (in some 
appropriate gauge), 

\im,,y^^^r, 9) = |$/| exp {lO) (2.3) 

where r^ = x^ + y^,iaii{9) = y/x. The cosmic string runs parallel to the z axis and is 
a static solution of the equations of motion. Apart from the cylindrical symmetry, the 
solution is invariant under translations in the z direction. 

The translationally invariant and homogeneous false vacuum at $ = $/ is metastable 
and its decay rate is determined by (|2.1| ). The corresponding bounce configuration $ is 
likely to be an 0(4) symmetric field configuration and the tunneling rate corresponds to 
the nucleation rate of unseeded bubbles. However the lowest energy state with an infinite 
cosmic string is also a local minimum of the energy functional and one can estimate the 



tunneling rate associated with its decay. We will treat a region of space with an isolated 
straight cosmic string as a classical false vacuum which is simply non-homogeneous. The 
formalism behind equation (|2.1| ) applies. The only important difference is that the field 
configuration at the false vacuum is a space dependent function which is not invariant 
under translations in x and y directions. The bounce configuration in this case is likely to 
be more complicated and will certainly lack the 0(4) symmetry of the usual bounces. The 
corresponding tunneling rate will be interpreted as the nucleation rate of critical bubbles 
along the length of the cosmic string. 

When the false vacuum is chosen to have a single infinite cosmic string with the above 
properties, the Euclidean action for the new bounces must be computed with the energy 
density of the false vacuum normalized to zero. Let us denote the field configuration 
corresponding to the local minimum of energy with a single cosmic string by $°, A^ 
and the corresponding covariant derivative and field strength by D^^^ and F° . The 
appropriately normalized Euclidean action is then (dropping the subscript on S): 

S = Ki + K2 + P + Fi + F2 + F3, (2.4) 

With/ri = /rf^xi [\D,^\^ + |A$|'] ; K2 = !d^x\ [|D,.$|2 + |Dj,$|2 - |D,,$0|2 _ |Dj,$0|2] ; 

P = !d'x [[/($) - f/(<|.°)] ; Fi = / dS^ fe - (F° )2l ■F2 = ! d'x^ \Fl + Fl + Fl + F, 



and F3 = Jrf^x^ Witl- Here we have used the result that, •l'^ and A^ have cylindrical 
symmetry around the z axis with A^ = A^^ = K^ = F°^ = F^y = F°, = F^^y = 0. We will 
refer to this solution of the equations of motion as the trivial solution G^. The regularized 
quantities K2, Fi and P are finite in the neighbourhood of G°. 

The search for stationary points of the action is greatly simplified by imposing sym- 
metries on the field configurations. We will impose rotational symmetries in the x, y as 
well as in the z,t plane. Introducing the polar coordinates E? = z"^ + t"^, tan(7) = z/t, 
r"^ = x^ + y'^ and tan(6') = y/x we write $(x, y, z, t) = a{r, R)exp{i9), Aq = v{r, R)/r and 
Ar = Ar = A^ = where a and v are real. The action reduces to: 

S = 2n^J JRdRrdr ^ (^\dRa\' + \dra\') + ^{1 - vf^ + j\^{drvf + U{a) - Co , 

(2.5) 
where the counterterm Co is obtained by replacing a and v by the two functions ao{r) and 
Vo{r) corresponding to the trivial solution G^. The above ansatz for the fields allows us 
to set Ar, Aft, A^ = 0. Variations in $ and A that do not preserve the ansatz contribute 
only positive definite terms to the action and therefore the ansatz minimizes the action 



with respect to those variations (this can also be seen as a consequence of the principle 
of symmetric criticality 0). Consequently a stationary point of the action in (|2.5| ) is a 
stationary point of the full action ( p.4|) . The ansatz also implies that ^2 = ^3 = 0. In 
the rest of the paper we always take At,Az,Ar,F2,Fs = and enforce the 0(2) x 0(2) 
symmetry described above. Also, by fixing the phase of $ in the above fashion we fix the 
gauge completely. 

The trivial solution G° has the same 0(2) x 0(2) symmetry and is obtained by finding 
the minimum energy solution with the boundary conditions (|2.3| ). The energy functional 
is given by 
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Minimizing E gives the equations of motion 
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G" is the solution ao(r), t>o(r) of the above equations that satisfies the boundary conditions: 
ao(0) = 0, ao(oo) = |<l>/|, t;o(0) = 0, t;o(oo) = 1. 



Note that with the normalization of (|2.4| ) the action of G° is zero and it is a local 



minimum of S. In contrast, a bounce solution in this case must have the following 

properties: 

(i) It starts at Euclidean time t = —00 from the static cosmic string solution and ends 

at Euclidean time t = +00 at the same solution but is not given everywhere by the the 

trivial time translation of the static string solution. 

(ii) The solution, which extremizes the action, is a saddle point of the action with a single 

direction of instability, i.e with a single negative eigenvalue of the Hessian grf-fr^- 

(iii) All velocities {d^/dt and dA^/dt) are zero for a temporal slice at some finite time 

(which we can choose to be at t = 0). The field configuration at t = is often called the 

turning point of the bounce. 

The direction of instability associated with the bounce configuration is intimately 
related to scale transformations of the fields. For a systematic exploration of these trans- 
formations it will be useful to define a one parameter deformation of the action (|2.4|) in 
the following manner (after dropping F2 and -F3): 

S[r,] = K,[r,] + K2 + P + Fr- (2.8) 



where K21 -P, Fi are as defined before, Ki[ri\ is a deformation of Ki defined by : 

K-i^[r]\ = fd^x- [\D,<^\^-" + \Dt<l>\'^-''^] . (2.9) 

In tlie limit rj ^ we recover tlie usual Abelian Higgs model (with the constraint 

A,,Ar,F2,F3 = 0). 

Under the scale transformation {z,t) -^ {Xz,Xt), with A being a positive number, 
the terms in the deformed action Sit]] transform as: Ki[ri] -^ X^Ki[ri] ■,{K2,P,Fi) — *• 
X'^{K2,P,Fi). Under these transformations the configuration G^ is left invariant while 
the finiteness of the action is preserved in the neighbourhood of G^. An extremum of 
the deformed action must be invariant under the scale transformation. Therefore at a 
stationary point of the deformed action we have, 

6S[v] 
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The second derivative of S[r]] at this point is: 



^_^ = 7] (KiM) +2[P + K2 + F^] = 0. (2.10) 

A — 1 
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7]{r] - l)K,[r]] +2[P + K2 + F^] . (2.11) 



Since Ki [77] > for any non-trivial field configuration, the above quantity is negative for a 
non-trivial extremum of 5* [77] if 2 > r] > 0. We will define the quantity V = [P + K2 + Fi] 
and prove the following statement. 

Statement 1: A bounce configuration exists in the Abelian Higgs model with the 
deformed action provided 2 > r] > 0. The bounce is the global minimum of the action in 
an open space of field configurations. 

Proof: Let the space of all smooth field configurations having the boundary condition 
given by property (i) of the bounce, satisfying ( p. 10 ) and having the 0(2) x 0(2) symmetry 



^{x,y,z,t) = \^\{r, R)exp{i6) be C. The trivial solution G^ lies in C Every point in G 
automatically satisfies properties (i) and (iii) of the bounce. We need to show that there 
is a point in G that also satisfies property (ii). 

We can split G into two parts G = Go (B Gi where Co is a connected space in G 
containing G°. In fact Gq consists of a single point G" since there exists a ball around G^ 
such that all field configurations contained in the ball except G° have V > 0. This follows 
from observing that at the false vacuum the energy density is minimized and therefore G^ 
is a local minimum of S[r]] with V = 0. However all nontrivial points on G have Ki[r]] > 
and therefore by ( p.lO|) V < 0. 
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The space Ci is not empty since one can always construct a non-trivial field config- 
uration that satisfies (|2.1CI|) and the correct boundary conditions. By construction, the 



global minimum G in Ci is is an extremum of the action with respect to all allowed vari- 
ations. By (|2.11|) it also has a single direction of instability corresponding to the scale 
transformations {z,t) -^ X{z,t). Therefore it has property (ii) of the bounce. Q.E.D. 

The limit rj = can be approached in the following way. By definition, G is the global 
minimum of Ci. First we show that G remains a maximum with respect to variations in a 
direction "orthogonal" to Ci even when rj = 0. The easiest way to see this is to consider a 
curve C\p] with p G [—1, 1] in the space of all field configurations with the right boundary 
conditions and passing through G (C[0] = G) in a direction "orthogonal" to the space Ci. 
For 2 > 1] > one can draw this curve so that one end of the curve is at G^ {C[—l] = G^) 
and the other end approaches the point G^ where $ = $t (^[1] -^ G^)- Note that the 
space Gi forms a boundary between the neighbourhoods of G^ and G^. This is because 
l^[G'°] = and ^[G^] -^ — oo, therefore the curve C[p] must intersect the codimension 
one surface Gi where V^ < 0. The curve can be drawn so that Spfp]] < S[G] for all 
p y^ 0. As r/ ^ 0, this property of the curve remains true and G remains a maximum of 
the action with respect to the variation that corresponds to moving along this curve. In 
this sense G remains a saddle point even when rj = 0. In principle the negative eigenvalue 
of the Hessian gf^ ^,2 may approach zero as r/ —>■ 0. Nevertheless, the zero eigenvalue 
can not be generic since there is no underlying symmetry that requires a zero eigenvalue. 
Phrased differently, all eigenvalues of the Hessian are likely to be of the order of m"^ where 
m is some characteristic mass scale in the theory. An eigenvalue can not be made zero 
without fine tuning the parameters of the theory, unless there is a symmetry to protect 
its smallness. Hence, for generic theories the limit rj ^ can be safely taken. 

Although the above discussion has been carried out in the context of an Abelian Higgs 
model with a single Higgs, everything said so far applies to the case of multiple Higgs 
fields. The generalization to the non-Abelian case is straightforward too; one simply 
needs to reduce the action in the non-Abelian case to an action resembling the Abelian 
model by dropping all components of the gauge field except one corresponding to a broken 
generator that generates a non-contractible loop on the vacuum manifold. 

So far we have proven the existence of a bounce with a large symmetry (0(2) x 0(2)) 
that that also happens to be the global minimum of the action in a simple space of field 
configurations. As we shall see, the latter property makes it quite simple to bound the 
action of this bounce from above. If this bounce is the bounce with the least action 



then its action determines the tunnehng rate of equation ( ^.ll ). If there is a bounce with 
a smaller action (and a lesser symmetry) then G provides only an upper bound on the 
least bounce action and consequently a lower bound on the tunneling rate. This bound 
itself may be useful in phenomenological applications and in section 4 we suggest simple 
numerical techniques for obtaining this bound. 

3 Phenomenological Applications 

The formal expression for the tunneling rate is similar to equation ( p.l| ). However we 
must drop the zero modes from the primed determinant. There are just two zero modes 
of the bounce solution coming from translations in the z and t directions. To leading 
order in ft we have: 



r (S[^\ /^-5[<l>]^^ det'[-d^ + U"{^,A)] '/' 



^ = L = U^r"Pl n 



det[-92 + f/"($o,AO)] 



(3.12) 



where S is the action of the string induced bounce [$, v4] (if and when it exists). Note 
that the configurations ^^ , A^ are space dependent. Corresponding to the two zero modes 

/ — \ 1/2 

there are only two factors of ( i^ 1 on the R.H.S 0. The ratio of the determinants in 
( ^.12| ) now has the dimensions of mass squared and the transition rate is to be interpreted 
as the probability of transition per unit time per unit length of the string. Collecting all 
non- computable parts into a single prefactor A we can rewrite the above equation as: 

7 = >lexp(-S[$, A]//i) . (3.13) 

To a good approximation A can be replaced by the characteristic squared mass scale of 
the theory and only the exponent needs to be determined carefully. 



Phenomenologically the tunneling rate obtained from ( 3.13|) is useful in two distinct 
situations. The first is the case when the theory has a viable false vacuum and one requires 
that the phase transition should never occur. The second case is a GUT model at high 
temperatures where one would like to calculate the high temperature vacuum decay rate 
and see if the string induced bubbles can percolate and complete the phase transition. 
We will briefiy discuss both cases. 

Consider the case of the viable false vacuum first. In this case, the false vacuum should 
have a lifetime longer than the age of the universe. A constraint is therefore obtained 
from the zero temperature tunneling rate. Multiplying the transition rate by the total 
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world sheet area of the cosmic strings, one gets the expected number of bubbles seeded 
by cosmic strings in a given space time. Most of the contribution to the string world 
sheet area comes from later times when the string network has entered a scale invariant 
distribution [§. During these times about 80% of the string length is contained in one 
long string per horizon. The total world sheet area available for bubble nucleation is 
therefore ~ H^ where H is the present horizon distance. The false vacuum would be 
regarded unstable if x ^ ^^ > 1- If the mass scale of the theory is greater than the 
electroweak scale then this condition translates into the following bound on the bounce 
action: 

5[$,A]<200^. (3.14) 

The same arguments repeated for unseeded bubbles yields the condition: 

5[$o,lo] <400/i. (3.15) 

Here [^^jA*^] is the ususal 0(4) invariant bounce and we have used the condition that 
^ X H^ > 1. Notice that the tunneling rate is expressed as the number of bubbles per unit 
time per unit volume and correspondingly one uses the total space time volume available 
for nucleating the bubbles. 

Cosmologically, there is no inflation in the false (but viable) vaccum and the bubbles, 
once formed, eventually percolate. The stability of the false vacuum therefore demands 
that both conditions ( |3.14| ) and ( p. 15] ) be untrue. The significance of this lies in the 
fact that a-priori the bounce actions S'[$, A] and 5'[$ ,A] are not related in any simple 
manner and constraints obtained on the parameter space of the model by ( |3.14D are 
independent and may be stronger than the constraints implied by ( p.l5| ). The simplest 
way to see this is to note that it is possible to have theories where the cosmic strings are 
classically unstable |]l3l . By a continuous deformation of the parameters of a theory one 



can go from a regime of classically unstable strings to metastable strings. Clearly there is 
a continuous family of theories in between where the bounce action in ( |3.14]) is arbitrarily 



small compared to the bounce action in (|3.15|) and the string induced tunneling is the 
dominant tunneling effect. 

Now consider the case when the viable vacuum is the true vacuum. The interesting 
case arises when the universe gets trapped in the false vacuum and supercools. Then 
the phase transition must be completed by seeded and unseeded bubbles in an inflating 
universe which (in comoving coordinates) is described by the Robertson- Walker metric: 



dr' = dt'-R\t) 



l-kr'^^''"''' 



(3.16) 



where k = —1,0 or +1 depending on whether the universe is open, flat or closed respec- 
tively. The scale factor R{t) grows according to the equation: 

■in k , , 

rP-^, (3.17) 
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where Mp = 1.2 x 10^^ GeV is the Planck mass. The energy density p is given by 
P = Po + (fo ) vT'^ where po is the energy density of the false vacuum, rj is the number of 
effectively massless degrees of freedom and T is the temperature. When the temperature 
is below the critical temperature Tc where the false vacuum becomes metastable, the 
energy density of the vacuum po niakes the dominant contribution to p. Approximating 
the energy density p by the vacuum energy density po we obtain in the period of inflation: 
"r ~ XPo; and the scale factor grows exponentially with time {R oc exp(xt)) where x = 

/ \ 1/2 

( llF ) ■ With adiabatic expansion we also have: RT = constant, so the temperature 
falls exponentially with time. 

In an inflating universe, the nucleation rate of unseeded bubbles grows as the volume 
of the universe (~ i?'^) while the same rate for seeded bubbles grows only as the length 
of the cosmic strings (~ R). Nevertheless, the seeded bubbles can play a significant and 
even the dominant role in the phase transition. To see this consider the probability that 
an arbitrary point remains in the supercooled false vacuum at time t. This was found in 
ref. to be: 





p{t) = exp 


in 
3 


I ft dt" 1 

[Jf R{t")\ 


3 

is the cc 



dt'n{t')V{t,t') 



(3.18) 



where V{t,t') = ^ J^., -^nrn. is the coordinate volume at time t of a bubble formed at 
time t' and n(t')dt' is the number of bubbles formed per unit coordinate volume between 
the times t' and t' + dt' and t^ is the time corresponding to the temperature Tc. If the 
tunneling rates associated with seeded and unseeded bubbles are respectively 70 and 7 
repectively, we can write: 

n{t) = 7oit)R'{t) + -i{t)aR{t)R\tc) , (3.19) 

where a is the total string length per unit coordinate volume at time t^. 

Both 7o and 7 are the finite temperature tunneling rates which depend on time (or 
temperature) and can be evaluated by computing the bounce action in the 3 dimensional 
Euclidean space using the finite temperature effective potential. The tunneling rates will 
however be greater or equal to their zero temperature values 7(0) and 7o(0). With the 
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approximations made above we can re- write equation ( |3.18D as a function of temperature 



p{T) < exp 



Tc dJ\ /7o(0) a7(0)\ / f^' 



(3.20) 



with d = j^. One can look at the above expression in two hmiting cases. If the seeded 
bubbles have little effect on the phase transition we can drop the terms multiplying 7 
from ( |3.20|) . Then in the limit T — i> we have: 



piT ^ 0) < exp 



-rf7o(0)log ^ ^ 



(3.21) 



If (i7o(0) is of the order of unity or greater, a fast phase transition is accomplished. 
For small values of (i7o(0) ((^70(0) < 10~^ may be sufficiently small) there may be no 
percolation or the resulting universe may be too inhomogeneous if reheating is achieved 
through bubble collisions |T^. Therefore viability of a model requires that either the 



transition be fast with only a mild inflation or the first order phase transition is followed 
by a slow roll-over of the infiaton as in new infiation. 

The effect of the seeded bubbles can be seen by putting 7(0) = in (|3.20|) . Then we 
get: 

p{T ^ 0) < exp [-rf(T7(0)] . (3.22) 

In this case the fraction of unconverted space approaches a non-zero value {p 7^ 0) as 
T ^ 0. However if (icr7(0) > 1 then the phase transition is fast and the resulting universe 
is homogeneous with the temperature at the end of the phase transition being of the 
order of Tc- Comparison of ( p.21|) and (|3.22|) shows that for (77(0) > 7o(0) the seeded 



bubbles make a significant contribution to the phase transition. The value of a depends 
on the initial distribution of the cosmic strings. There are two extreme cases. If the 
temperature at which the strings are formed is close to Tc then the string network is at 
the friction dominated period when the infiation begins and a r^ T^. If on the other hand 
the strings are formed at a much higher temperature than T^ then they would have a 
scaling distribution at T^ and we would expect about one long string per horizon. This 
implies that a ~ T^/M^i. The characteristic mass scale associated with the tunneling 
rates 70 and 7 is likely to be Tc. Therefore we can write: 70 (0) ~ T^exp{—S'^/h) and 
7(0) ~ T^exp{—S/h), where S and 5* are the corresponding zero temperature bounce 
actions. Using these expressions the condition for the seeded bubbles to be significant in 
the transition may be stated as: S — a < S^, where a lies in the interval [0, log(T^^/Mp;)]. 
Given the fact that the two bounce actions are independent of each other, in any realistic 
model the above condition is as likely to be satisfied as not. 
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4 Computing The Bounce Action 

In this section we present a simple technique for actually computing the the zero tem- 
perature tunneling rate. The basic idea is similar to the suggestions of ref. ^ and ref. 
W^ . Instead of looking for a saddle point of the action in the space of all configurations, 
one looks for the global minimum of the action in the reduced space Ci. The restriction 
to the reduced space can be done by Lagrange multipliers. We may define an improved 
action 5* by: 

S = S + EiMi . (4.23) 

The extra terms Mi added to the action are of the form CjE"" where Cj is a Lagrange 
multiplier, n is a positive number and Bj are functionals of the fields that vanish on the 
reduced space. For instance, S = S + CiV^ has local minima at the trivial solution G° as 
well as G. By choosing ci large one can essentially force an iterative minimization to take 
place along a path on the surface Ci. Clearly the minimum one may obtain is either G or 
some other local minimum G' of S* in Ci. In either case one reaches a bounce. However, 
since Gi is an infinite dimensional space, in practice one can not arrive at the point G or 
even a local minium G' through finite methods. Indeed, the result of the search is likely to 
be a point which lies close to Gi but may not even be on Gi. Nevertheless, by performing 
a small scale transformation {z,t) -^ X{z,t), one can make the final point of the search 
lie exactly on Ci. Thus irrespective of the accuracy of the numerical method, the search 
always succeeds in finding an upper bound on the action of the bounce G which sets a 
useful lower bound on the tunneling rate. 

Another advantage of this method is that it admits extreme simplifications. Since 
one can at best hope to find an upper bound on the action of G, the computation can 
be simplified enormously at an early stage by making the search space small. The only 
required element of the computation is that condition V^ = be enforced so that the result 
of the search is a point on Gi. We will use the example of the Abelian Higgs model to 
illustrate a simplified search of this kind where the problem will be reduced to solving an 
ordinary differential equation in a single variable. (As mentioned before, the non-Abelian 
case is reducible to the Abelian one). We start with the action given by ( p.5| ). The most 
interesting case arises when the potential f/($) has a local minimum at |<l>| = $/ 7^ 
and a global minimum at $ = 0. This is, for instance, the case in the MSSM if the true 
vacuum is viable and the false vacuum breaks hypercharge. The cosmic strings in the 
false vacuum may be metastable or even classically unstable with respect to an expansion 
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of their cores. In the true vacuum, they simply do not exist. 

The first step is to determine the functions ao(^) and fo(^) numerically by minimizing 
the energy functional ( p.6| ) or by solving the equations ( |2.8| ) numerically. Next, consider 
the field configuration a{r, R) = ao{r — r'{R)) for r > r'{R) ; a{r, R) = for r < r'{R) ; 
v{r,R) = t'o(r — r'{R)) for r > r'{R) ; v{r,R) = for r < r'{R) ; where r'{R) is a real 
positive function of R satisfying r'{R)\ii^oo = and ^|i?=o = 0. The field configuration 
consists of a string that nucleates a bubble of true vacuum in its "core" at some t < 0. 
The bubble expands as one approaches t = and then it contracts and vanishes at some 
t > 0. For a given value of z and t, the cross section of the bubble in the x, y plane is a 
disc of radius r'. 



1 

r'(R) 



2r- 

\ / 

V 



R- 



FIG. 1 



FIG. 2 



FIG. 1. A typical form for the function r'{R). The function goes to zero as R 



CXD. 



FIG. 2. A schematic diagram for the string "core" after the nucleation of the 
bubble. The radius of the bubble in a plane perpendicular to the z axis 



IS r 



The above construction gives us a space of field configurations satisfying the correct 
boundary conditions and depending on a single real function r' of R. Restricted to this 
space of field configurations, the action is simply a functional of r'{R) which can be 
thought of as the Lagrangian of a particle with coordinate r' integrated over time R. We 
have: 



S 



2n^ / RdR 



(Li + r'L2)idnr'f + {L, + r'U) + ^r'^ 



Cn 



(4.24) 



where L 



^^{l-v,f + ^M'Vo?]: L, 



Jrdr{draoY] L2 = J dr^draof; L3 
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/ r drL\ and L4 = / drL. The integrals in r must be done numerically with the known 
functions Oq and Vq. Then the bounce in the reduced theory is obtained by solving the 
equation of motion implied by ( [4 .241 ) which is an ordinary differential equation. This is 



exactly the same computation that one needs to do to evaluate the tunneling rate in the 
usual translationally invariant case for a theory with a single real scalar field. It can be 
easily done by the shooting method that has been described elsewhere 0, P, [l^. 



Note that at the extremum, the action is extremized with respect to scale transforma- 
tions of R which correspond to scale transformations in the z, t plane. Thus the extremum 
automatically lies on the surface Ci and gives an upper bound on the action of G. In this 
case there is no further need to include Lagrange multiplier terms as in (^4.23 ) to force the 



result to lie on the surface Ci . This feature is typical of reductions of this kind where the 
final action depends on a single real function 0. 

5 Conclusions 

We have considered vacuum tunneling from non-homogeneous false vacua. A summary 
of our results is as follows: 

(i) There are new tunneling phenomena associated with topological defects in a false 
vacuum. The rate of vacuum tunneling through bubbles seeded by defects is independent 
of (and can be arbitrarily large compared to) the tunneling rate through bubbles that are 
not seeded by the defects. 

(ii) We have identified a bounce configuration (and call it G) which allows a semi- 
classical estimation of the tunneling rate. 

(iii) It is possible to search for G numerically, by using simple minimization of a 
function on the lattice. We suggest a method of reduction by which one can do the search 
over a reduced space and find an upper bound on the action at G. The corresponding 
tunneling rate is a lower bound on the actual tunneling rate, and may be useful in placing 
constraints on the parameter space of realistic field theory models. 
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